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Abstract 

Cn . It is known that the moments of the maximum value of a one-dimensional condi- 

tional Brownian motion, the three-dimensional Bessel bridge with duration 1 started 
from the origin, are expressed using the Riemann zeta function. We consider a system 
of two Bessel bridges, in which noncolliding condition is imposed. We show that the 
moments of the maximum value is then expressed using the double Dirichlet series, or 

'^ I using the integrals of products of the Jacobi theta functions and its derivatives. Since 

^2 ' the present system will be provided as a diffusion scaling limit of a version of vicious 

walker model, the ensemble of 2-watermelons with a wall, the dominant terms in long- 
time asymptotics of moments of height of 2-watermelons are completely determined. 

CN ' For the height of 2-watermelons with a wall, the average value was recently studied by 

J^ . Pulmek by a method of enumerative combinatorics. 

C^ ■ Keywords Bessel process • Bessel bridge • noncolliding diffusion process • Riemann 

zeta function • Jacobi theta function • double Dirichlet series • Dyck path • vicious walk 

O ■ 1 Introduction 



Let B{t) = {Bi{t), B2(t), B3{t)),t > be the three-dimensional Brownian motion (BM), in 
which three components Bj{t),j = 1, 2, 3 are given by independent one-dimensional standard 
5^ i BMs. The three-dimensional Bessel process (BES3), X{t), started from x > is defined as 
the radial part of B{t), 

Xit) ^ \B{t}\ 

= Vi?i(t)2 + 52(t)2 + 53(t)2, t>0 
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Figure 1: Sample path of three-dimensional Bessel bridge with duration 1. 



with X{0) = X. BES3 is a diffusion process on R_|_ = {x G R : x > 0}, where R denotes 
the set of all real numbers. By Ito's formula we can show that it satisfies the stochastic 
differential equation of the form 

dX{t) = dB{t) + -TFTrdt, t > 0, X(0) 



X{t) 



X, 



where B{t) is the one-dimensional standard BM different from Bj{t)^s used to give B{t) 
above. We can prove that X{t) ^ cxd in t — > 00 with probability one for all a; > 0, i.e. BES3 
is transient. For the basic properties of BES3, see, for example, 3.3 C in [13], VI. 3 in [23] , 
IV.34 in [6]. 

The three-dimensional Bessel bridge with duration 1 started from the origin, X{t),t G 
[0, 1], is then defined as the BES3 conditioned 



X 



X(0) = and X(l) 



0. 



Figure 1 illustrates a sample path oi X{t) on the spatio-temporal plane (t, x) G [0, 1] xR+. In 
[1] , a variety of probability laws associated with conditional Brownian motions are discussed, 
which are related to the Jacobi theta function and the Riemann zeta function. One of them 
is the probability law of the maximum value of X (t) ; 



Hi = max X{t). 



0<t<l 



Let E[i7^] be the s-th moment of Hi. The following equality is discussed in jl], 

s/2 



nHi] = 2[-) as) 



sgC, 



(1.1) 



;i.2) 



where C denotes the set of all complex numbers, and 



with the gamma function 



m = Isis - l)n-^/'T{s/2)as) 



T{s)= duu'-'e"', 3fJs>0, (1.3) 

'o 



and with the Riemann zeta function 

oo _ 
n=l 

See also Chapter 11 in [25] . 

We know the two facts; (i) the BM can be realized as the diffusion scaling limit of the 
simple random walk, (ii) the probability law of BES3 is equal to that of the BM conditioned 
to stay positive. Combination of them will lead to the following. For a fixed n > 1, consider 
one-dimensional simple random walks started from the origin, which visit only positive sites 
{1,2,3,---} up to time 2n and return to the origin at time 2n. Sample paths of such 
conditional random walks are called Dyck paths of length n in combinatorics. The height of 
Dyck path hi (2n) is defined as the maximum site visited by the walker. Let ( ■ ) denote the 
average over all Dyck paths with uniform weight. Then we will have the relation 

lto((M') = EM. ,.C. (1.4) 



The classical work of de Bruijn, Knuth and Rice in enumerative combinatorics [7] gives 

3 

(/ii(2n)) ~ i/vm — - + 0(1) in n ^ 00. (1.5) 



Here we should note that, through the relations (II. 2p and (II. 4p . if we only consider the 
dominant term in (II. 5p proportional to ^/n, this result in combinatorics means nothing but 
the fact ^(1) = 1/2. It is rather obvious if we know the following integral representation of 
^(s) due to Riemann, 

1 1 r°° 

as) = 2 + 4^(^ - 1) y du iu-'/'-' + u^'-^^/'-')i^iu) - 1), (1.6) 

where '(^(m) is a version of the Jacobi theta function 

00 
^(m) = J2 e~""'", u>0. (1.7) 

n=— 00 

Recently Fulmek reported a generalization of the result of de Bruijn, Knuth and Rice, 
by calculating the asymptotics of the average height of 2-watermelons with a wall [H] . In 



general, the uniform ensemble of A^-watermelons, A^ > 2, is a version of vicious walker model 
of Fisher [TU]. In this version the starting points and the ending points of A^ vicious walkers 
[i.e. nonintersecting random walks) are fixed to the sites located near to the origin. When 
we impose the condition to stay positive for all vicious walkers, we say "with a wall" (at the 
origin) [3l El EH [121 EDI [22] . The height of A^- watermelon is the maximum site visited by 
the vicious walker, who walks the farthest path from the origin. Let /i2(2n) be the height of 
2-watermelon with a wall. Fulmek showed 

3 

(/i2(2n)) ~ caV^ ho(l) in n ^ oo with Ca = 2.57758 ■■• . (1.8) 

Here the factor C2 = 2.57758 ■ ■ ■ of the dominant term proportional to ^/n was given by 
numerical evaluation of the "constant terms" in Laurent expansions of a version of double 
Dirichlet series. The terms are represented by integrals of functions expressed using the 
Jacobi theta function (ll.7p and its derivatives. It should be emphasized the fact that Fulmek 
succeeded in proving the N = 2 case of the conjecture of Bonichon and Mosbah [5], 



{hN{2n)) ~ v^(1.67iV - 0.06)2n in n ^ oo 

obtained by computer simulations for the average height hN{2n) of general A^- watermelons 
with a wall, A^ > 1. It seems to be highly nontrivial to extend his method to evaluate the 
asymptotics of higher moments {hi\i{2n)^), s > 2 ioi N = 2 and A^ > 3. See the paper by 
Feierl on the recent progress in this combinatorial method |9]. 

Here we propose a different method to calculate the dominant terms of all moments 
of height for 2- watermelons with a wall. We will perform the diffusion scaling limit first. 
Following the argument of fTEi [TU| [T2j , we can prove that the diffusion scaling limit of the 
A^- watermelons with a wall provides the noncolliding system of A^ Bessel bridges, X = 
(Xi(t), ■ ■ • , XN{t)) e W^ = {(xi, ■ ■ ■ , xn) : < xi < • • • < xn}, < t < 1. It implies 

,. //M2B)yv g(^,,_ ^^2. (1.9) 



n.^oo \ \^ ^2n 

where 

Hn= max XN(t). (1.10) 

0<t<l 

In the present paper we determine E[if|] for arbitrary s for the two Bessel bridges with non- 
colliding condition. Noncolliding diffusion particle systems are interesting and important 
statistical-mechanical processes, since they are related to the group representation-theory, 
the random matrix theory, and the exactly solved nonequilibrium statistical-mechanical mod- 
els {e.g., ASEP and polynuclear growth models) [19j. The present system of noncolliding 
Bessel bridges is related to the class C ensemble of random matrices discussed by Altland and 
Zirnbauer P,[2] (see Sect.V.C of [E]) and it is a special case with parameters (z/, k) = (1/2, 3) 
of the noncolliding generalized meanders [18] (see also [21] ) . 

As demonstrated in [TU], the noncolliding diffusion processes can be regarded as the 
multivariate extensions of Bessel processes. The present study suggests the possibility that 



the connection between the conditional BMs and the number theoretical functions {e.g., the 
Jacobi theta function, the Riemann zeta functions, and Dirichlet series) reported in [2^ H] 
will be extended to many particle and multivariate systems. 

The paper is organized as follows. In Sect. 2 we give the precise description of the systems 
and main results. In Sect. 3, we give proofs of our theorems for the N = 2 case and show 
formulas, which are useful to perform the numerical evaluation of moments. 

2 Models and Results 

2.1 Reflection Principle and Karlin-McGregor Formula 

The transition probability density of the one-dimensional standard BM is given by the heat- 
kernel 

1 r (y-xy 



pit, y\x) = , exp < > , X, -u G R, t > 0. 

By the reflection principle of BM, the transition probability density of the BM with an 
absorbing wall at the origin is given by 

pi{t,y\x) = p{t,y\x)-p{t,y\-x) 

v/2^V ' '^ ^ - 

If we put two absorbing walls at the origin and at x = /i > 0, then repeated application of 
the reflection principle determines the transition probability density of the absorbing BM in 
the interval (0, K) as 



oo 

v\{^,y\^) = ^ ^p{t,y\x + 2hn) -p{t,y\-x + 2hn)j 



n=—oo 



'2TXt 

n=— OO 



exp <j -l^Xv - (a; + 2hn)f I - exp <^ ~2t*^^ ~ ^~^ ^ 2hn)f 



for x,y & {0,h),t > 0. Since BES3, X{t), is equivalent with the BM conditioned to stay 
positive, and this process is realized as an /i-transform of the absorbing BM with a wall at 
the origin (see, for example, [19]), we will see that 

P{H^<h)= hm 4^ (2.1) 

for (11. ip . The limit of (12.11) can be readily performed and we have 

00 
P{Hi <h)= Y, £-=^^'"'(1 - 4h^n^) 
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Figure 2: Sample path of two Bessel bridges with duration 1 conditioned never to colhde. 

and the probabihty density is obtained as 

d 



qiih) 



dh 



P{Hi < h) 



n=l 



-2h^n- 



2»,2 



{4h'n^-3hn'). 



The s-th moment of Hi is defined by 

E[Ht] 



dhh''qi{h), 



and (ll.2p is derived, for which the following equalities are useful, 



dhh^e 



s„-2h''n^ 



T{s + 1) 



2-(.+3)/2^-{.+i)r^^^ + l)/2), 
sT(s), ^s > 0. 



(2.2) 



For A^ = 2, 3, ■ ■ ■, the noncoUiding A^-particle system of Bessel bridges with duration 1, 
with all particles started from 0, is denoted by Xj^{t) = (Xi(t), ■ ■ ■ ,Xj^{t)), where 

0<Xi(t) <X2(t) < ■■■ <Xiv(t), 0<t<l. 



The stochastic variable if at is defined as the maximum value of the iV-th Bessel bridge (11.101) . 
See Figure 2 for the N = 2 case. By the Karlin-McGregor formula [H], we will have 



P{Hn < h) 



lim Fh{yi,y2, 

Xj^O,yj^O,l<j<N 



,yN\Xl,X2, 



,Xn) 



where 

^<de^^[p2(l,%|a;fc)] 

Fh{yi,y2, ■■■, VnIxi, X2,---, xn) = ~^' ~ r i^ , — 77 (2.3) 

, detjpi (1,2/^1 Xfe)] 

for a; = (xi, ■ ■ ■ , x^\ y = (yi, ■ ■ ■ , y^) G Wj^ = {0 < Xi < • ■ ■ < xat < h}. The s-th moment 
of H]\f is then given by 

/•OO 7 

E[H-]^]= dhh'qNih) with qj^{h) = —P{Hn < h). 

Jo d,n 

2.2 Results 

Here we show our expressions for the moments of H2 of the two Bessel bridges conditioned 
never to coUide in t G (0, 1). First expression is given using the double Dirichlet series of the 
form 

(ni,n2)eZ2\{(0,0)} ^ 1 2; 

where Z denotes the set of all integers. 

Proposition 2.1 Let 

Z^(h) = T{a + 2b)Z{2b, 26; a + 2b). (2.5) 

Then 

nH^] = ^s[(s-l)(s2-2s + 12)Z,/2(0)-4(s + 4)(s + 6)Z,/2(l) + 64Z,/2(2) 

(2.6) 

Remark 1. As we show in Sect. 3.1, Zs/2{b),b = 0,1,2, have simple poles at s = and 
s = 2. By the prefactor s and by cancellation of the s = 2 poles among the three terms, 
however, the expression (12. 6p has finite limits in s -^ and s ^ 2. See Eg. (13. lip below. 



We can rewrite this result using the Jacobi theta function (11.71) and its derivatives, 
'&'{u) = d^{u)/du,r{u) = dH{u)/du^. 



Theorem 2.2 Let 



/OO 
duu''''-\d{uf-l}, (2.7) 



and 



Us) = ~^Us + A){s + Q) r duu''^+H\uf 



+ ((2 -s)+ 4)((2 - s) + 6) / duu^^-'^/^+^d'{u) 



Q />00 1 

+ \ / du (m^/2+3 + «(2-)/2+3)^//(^)2 ^ _L^(^ _ 2)^(1)2^ 

3 Ji 12 



:2.8) 



Table 1: Numerical values of moments 



s 





1 


2 


3 


4 


5 


EH^ 


1.0 


1.253314 


1.644934 


2.259832 


3.246969 


4.873485 


EH^ 


1.0 


1.822625 


3.395156 


6.463823 


12.576665 


25.005999 



Then 



E\m 



2 






se C. 



(2.9) 



Remark 2. By the integral representation (11.61) . it is clear that ^(s) satisfies the functional 
equation 

e(i-s) = e(s), SEC. 

It is interesting to see that the function ^2(5), which appears in the expression (12.90 . satisfies 
the functional equation 

W2-s)=Us), sec. 



As will be explicitly given in Sect. 3. 3, ^(s), Kq{x) and ^2(5) are expressed using series of 
the incomplete gamma functions. Numerical evaluation of the incomplete gamma functions 
is easy, and the series converge rapidly. Actually we have readily obtained the values of 
moments for A^ = 1 and A^ = 2 as shown in Table 1. (The trivial result E[l] = 1 is 
obtained by setting s = in (II. 2p with (11.61) and in (12. 9p . Since we know Euler's work 
on the relation between ({2n),n = 1,2,3, ■■■ and the Bernoulli numbers, Eq. (ll.2l) gives 
E[H^] = C(2) = 7rV6 = 1.644934- ■ ■ and E[Hf] = 3C(4) = 7r730 = 3.246969 ■ ■ -.) By the 
relations (11.41) and (II. 9p with N = 2, from the values in the s = 1 column in the Table 1, 
the dominant terms of the previous results (II. 5p and (II. 8p are reproduced; 

{hii2n)) ^ 
{h2i2n)) ^ 



2n X E[Hi] 
2^x 1.253314' 
2^ X E[H2] 
2^x 1.822625' 



vrn, 



2.57758- --x ^/^. 



3 Proofs and Numerical Calculations 

3.1 Proofs of Theorems 

In this subsection we will prove Proposition 12.11 and Theorem I2.2[ Let rd{yi,y2\xi,X2) and 
rn{yi,y2\xi,X2) be the denominator and the numerator of the RHS of f l2.3p . respectively, 
given by 2 X 2 determinants, when N = 2. We have obtained the estimations 

rd(a;i,a;2|xi,a;2) = --xjxKxi - X2f{xi + X2f + 0{x\°,xl^), 

671 



oo oo 



rn{Xl,X2\Xi,X2) = —xlxl{Xi-X2f{Xi+X2y ^ ^ 6 2'^'('^i+"2)Q^(^^^ ^2) 

ni=— 00 n2=— 00 

+0{xl,xl), 

for xi, a;2 -C 1 with 

Qh{ni, ria) = 3 - ASh^nl + 72h^nt + 72h'^njnl - 32h%l - 96h%'lnl 
+128h^n^^nl - 128/i^n^n^. 

Then the following results are concluded. 

Lemma 3.1 

P{H2 <h)= Yl e-2^'("?+"') A;,(ni, ns) 



2 



with 



And then 



(ni,n2)eZ 



32 

Ah{ni, na) = 1 - Wh^nl + 24/i^n^ + 2Ah'^njnl h^n^ - 62h^n{nl 

3 ^ ^ 3 ^ ^ 

(ni,n2)GZ2\{(0,0)} 



o , 

+16/i^n^ + lQ2h^nlnl - mh^n\nl - QAh^nlnl + 64/z*n5n^|. 

Proof of Proposition \KT\ In order to describe the moments E[/72] = /o dhh^q2{h), we 
introduce the following notation, 

rf/i/i^+^-i+^e-2'^' ("?+"'). (3.1) 

ini,n2jfcZ,-\UU,u;| 



Then, by Lemma [3 -H we have 



E[K 



15/(2, 0) + 60/(4, 0) + 60/(2, 2) - 60/(6, 0) - 180/(4, 2] 
+16/(8, 0) + 192/(6, 2) - 80/(4, 4) - 64/(8, 2) + 64/(6, 4)|. (3.2) 



Next we rewrite /s(a, (3) using the Gamma function 01.31) . In the integrals in (13. ip . we change 
the integral variable from h to u hj u = 2h'^{n\ + n^), respectively. Then we have 

/,(«, (3) = 2-("+'^+2+s)/2p^(^ ^p^ s)/2)Z{a, (3- {a + (3 + s)/2), 
where Z(a,/3;7) is defined by (12.41) . From (13. 2p we will see that 



nm 



j^2{2-s)/2 

3 



151(1 + s/2)Z(2,0;l + s/2) 

+30r(2 + s/2){Z(4, 0; 2 + s/2) + Z(2, 2; 2 + s/2)} 

-151(3 + s/2){Z(6, 0; 3 + s/2) + 3Z(4, 2; 3 + s/2)} 

+2r(4 + s/2){Z(8, 0; 4 + s/2) + 12Z(6, 2; 4 + s/2) - 5Z(4, 4; 4 + s/2)} 

-4r(5 + s/2){Z(8, 2; 5 + s/2) - Z(6, 4; 5 + s/2)} . (3.3) 

By definition ([231), 

Z(2.0;l + ./2) = l ^ ^i±!i_= lz(0.0;,/2), 

(ni,n2)eZ2\{(0,0)} ^ 1 27 



Z(4,0;2 + s/2) + Z(2,2;2 + s/2) 



= -Z(0,0;s/2) 



n^ + 2nln\ + ra| 

(ni,n2)eZ2\{(0,0)} ^ 1 27 



E 



Z(6, 0; 3 + s/2) + 3Z(4, 2; 3 + s/2) = - 



-Z(0,0;s/2), 



(ni,n2)GZ2\{(0,0)} 



nv 



3n'in2 + 3nin: 



2^4 
2 



nn 



(n? + n2)3+«/2 



Z(8, 0; 4 + s/2) + 12Z(6, 2; 4 + s/2) - 5Z(4, 4; 4 + s/2) 

Til + 12n\n\ — bn\n\ 



E 



(ni,n2)eZ2\{(0,0)} 

i E 

2 ^^ 

(ni,n2)eZ2\{(0,0)} 



(^2 + n2)4+«/2 

(n? + n2)4 + Mlnl{nl + n^)^ - 32n\n\ 
(nf + ni)4+«/2 



-Z(0, 0; s/2) + 4Z(2, 2; 2 + s/2) - 16Z(4, 4; 4 + s/2). 



10 



n^rio — n^nk 



Z(8.2i6 + ./2)-Z(6,4;5 + ./2)= ^ JlJ^Z 

(ni,n2)GZ2\{(0,0)} ^ ^ '^' 

1 ^^ n\n\{n\ + n^Y — An\n2{n\ + nl) 

~ 9 2-^ 



(ni,n2)eZ2\{(0,0)} ^ 1 2; 

= ^Z(2,2;2 + s/2)-2Z(4,4;4 + s/2). 
Using the above results, fl3.3p is rewritten as 

E[/72l = ^-^^ [ci(s)Z(0, 0; s/2) + C2(s)Z(2, 2; 2 + s/2) + C3(s)Z(4, 4; 4 + s/2) 

with 

ci(s) = -yr(l + s/2) + 15r(2 + s/2)-yr(3 + s/2) + r(4 + s/2) 

= l^s{s-l){s^-2s + l2)T{s/2), 
lo 

C2(s) = 8r(4 + s/2)-2r(5 + s/2) 
= -sr(s/2 + 4) 

= -is(s + 4)(s + 6)r(s/2 + 2), 

C3(s) = -32r(4 + s/2) + 8r(5 + s/2) 
= 4sr(s/2 + 4), 

where (12.21) has been used. Then Proposition 12.11 was proved, i 

Let 1{^} be the indicator function of condition cu; 1{^} = 1, if cj is satisfied and Ij^^} = 0, 
otherwise. The following equality is derived. 

Lemma 3.2 

Za{h) = vr" P du u^+^'-' \(J-{}{u)\ - l|,=o} 1 • (3.4) 

Proof. By definition (12.51) 

Za{h)= E ^,2,\L^2. r(a + 26). (3.5) 

(ni,n2)6Z2\{(0,0)} ^ ^ '^' 

By changing the integral variable m by w with u = vr(ra^ + 722)10 in the integral (II. 3p . we have 

r(s) = TT'inl + nlY / dw u;^-ie~"("?+"2)-. 



11 



Then (13.51) becomes 

Za{b) = TT 



"+2^ ' dww 





a+2b-l 



E 

(ni,n2)GZ2\{(0,0)} 



n^^n^^e"''^"i+"2)«' 



TT 



a+2b 



POO I / °° 



n2b^-nn\ 



L{fe=0} 



Since (—Trn ) e~ 



dw 



-e"™''", (IS3D is obtained. 



Then we have the following expressions of moments 
Proposition 3.3 

1 /7r\s/2 

^ 2J 24 V 2 



{s - l)(s' -2s + 12) / duu'/''-^{^{u)^ - 1} 

Jo 

/'CO /"CXD 

-4(s + 4)(s + 6) / duu''^+^d'{uf + QA duu''^+^d"{uf 
Jo Jo 



(3.6) 



Proof of Theorem \2.S[ By the reciprocity law of the Jacobi theta function [1] 



^(u) 



y," 



u \u 



, 3fJM > 0, 



we can show the following, 



h 



duu''^~\^{uf-l} 





r) r) /'OO /'OO 

-- + ^-^ + / duU-''^{^{uf - 1} + / rfMM^/2-l|^(^)2 _ ^1^ 



S s-2 



(3.7) 



(3.8) 



rfMM^/2+l^//^N 



-+ / duu'/^+^^'{uf+ / duu^'/^+^^{u)^'{u) 

2(s — 2) Ji Ji 

/CO -I /'OO 

du u-''^+H\uf + ^ / rfw m-^/'{^(m)2 - 1}, 



(3.9) 



h = I duu'/^+^"{u) 

'0 



12 



/CO /"OO 



{s-2) 



1 

oo 



+? r duu-''^+^d{u)d'{u) + Q I duu-''^+^d\u) 



9 



+— / duu-''^{d{uf-l]. 
16 Ji 



(3.10) 



The derivations are given in Appendix A. They show that /i has only two simple poles at 
s = (with residue —2) and at s = 2 (with residue 2), I2 does only one simple pole at s = 2 
(with residue 1/2), and /s does only one simple pole at s = 2 (with residue 9/8). Using 
them in the representation 03.61) . we have 



E[i7| 



1 /7r\^/2 
24 V2 



2(s^-14s + 12) 



+s|(s - l)(s' - 2s + 12)Kq{s) - 4(s + 4)(s + 6)7^1(5) + QAK^^s) 
+s{s - 2fKo{2 - s) - 4(s^ + 10s - 120)Ki(2 - s) + UK2{2 - s) 

-4(s^ + 10s-48)Ji(s) + 96J2(s) + 384J3(s)| , (3.11) 



where 



duu'''''-^{^{uf -I}, 

X) 

duu''''''H\u)\ 

X) 

dnn^/2+¥"(M)^ 

1 

XD 
XD 

duu"'''^''d{u)d"{u) 

X) 

rfMM-^/2+3^'(M)r(M). 



i^o(s) 

K2{S) 

Ms) 

J2is) 

Ms) 



Note that the singularities in /i , /2 , I3 have been all cancelled and all integrals converge for 
all s G C. 

Now we perform partial integrations to evaluate Ji(s), ^2(5) and J3(s). For Ji{s), 



Ms) 



u 



-s/2+1 



TOO /""O d / , \ 

d{u)^^ - / rfu — (^m-^/2+^^(m)J^(u) 



-l 00 1 

u-''^+'{}{uf + (s-2)Koi2-s) 

J 1 .^ 



M 



-s/2+1 



Jl(^). 
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Since 
we have 

For J2(s), we see 

J2(S) 



U-^/2+l^(^)2 



U 



-s/2+1 



J 1 



-Hif-1}, 



.Ms) = \is-2)Ko{2-s)-^{^iir-l}. 



(3.12) 



du 



y [ du^\u-i--H[u)]4,m 



du 



du 



-^(l)^'(l) + -{s - 4) Ji(s) - Ki{2 - s). 



Inserting (13.121) . we have 
Us) = -d{l)d'{l) --^{s- A){d{lf - 1} + l(s - 2)(s - A)K,{2 -s)-K,{2- s). (3.13) 



Similarly by partial integration, we obtain 



Ms) 



du 



+ -{s-6)K,{2-s)-Ms) 



= -l^'il)' + \is-6)K,{2-s). 
Using (IXT^ - dXTlD in (IHTTD . (1^:^ with ([ZTD and dZHD is obtained. 



(3.14) 



3.2 Expressions by Incomplete Gamma Functions 

Let r{z,p) be the incomplete gamma function of the second kind defined as 



T{z,p) 



duu^ e 



z~l „—u 



T{z)- f duu'-^e-'', ^z>0,p>0. 
Jo 



(3.15) 



As demonstrated in Appendix B, the integrals appearing in (II. 6p and in our result given by 
Theorem 12.21 are expressed using r{z,p). We have obtained the following expressions. 



E[Ht] 



2 



l + s(s-l) 



X 



7r-^/2 Y^ n-'r{s/2, Tin^) + n'/^-^/^ ^ n'-'r{-s/2 + 1/2, 



rm 



n=l 



n=l 



, (3.16) 
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and 



nm 



7r\s/2 
2 



l(l_,)(,2_2s + 12) 



X <^ 1 - 2sn-'^^ 



(oo oo _ oo , 

E E KT;i|jvir(.V2,x(„; + „!)) + J: -r(./2. 
ni=ln2=l ^ i 2^ n=l 



vm 



-4s{^(l)/(l) + 2(^'(l))2} + s6(s) 



(3.17) 



with 

Us) 



i{n-ns+^){s+6)j2Yl 



nfnl 



^ ^^ (n? + n2)^/2+2 

i=ln2=l ^1 2; 



r(s/2 + 2,7rK + n^)) 



oo oo 



2^2 



nrn: 



l'^2 



+7r^/-^((2 - .) + 4)((2 - 3) + 6) 5: 5: ^^ ^ /,^3 

ni=l 712=1 ^ 1 2/ 



r(-s/2 + 3,7rK + r2^)) 



32 



+T --'^EE 



i4^4 



i=ln2=l ^1 2; 



r(s/2 + 4,7rK + n2)) 



4 4 

nin2 



Z^ Z^ (^2^^2)-./2+5 
ni=ln2=l ^1 27 



r(-s/2 + 5,7rK + n^)) 



1 



+ __,(,_ 2)^(1)2. 



(3.18) 



The values in Table 1 were obtained by numerically calculating the incomplete gamma 
functions and the series of them with appropriate coefficients in the above expressions. 



Appendices 
A Calculation of /j, j = 1, 2, 3 

By the reciprocity law of the Jacobi theta function 03. 7p . 

/»1 /"OO 



duu^'^-^ 



1 



u 



+ duu'/^-^l i^i- 



U(-)) -A+ r duu'/''-\d{uf - 1}. 
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By calculating the first integral and by changing the integral variable u hj w = 1/u in the 
second integral, we obtain (13.81) . 
Set 

In the first integrand, we use the reciprocity law fl3.7p as 
d „, .\ Id /I „ /I 



du J ] du \ y u \u^ 

2 \uj \ udu \u 



4 VV^// \ujdu\uj \du \u 

By inserting this into the first integral in I2 and by changing the integral variables u by 
w = 1/u, we have 

/2 = - / dww'''^d{wf+ / dwW-'l^+^diw)-^diw) 

4 Ji Ji dw 



00 



2 



fX ^ / \ \ / - r. /o I 1 / (-*' 



+ / dww-''^^^ ( — ^(u;) J + / dww''^+^ ( —^{w) 



The first integral is rewritten as 



- / dww-''^d{wf = - dww-''^ + - / rfww^-^/'I^H' - 1} 
+ - / dww~'/^{d{wf-l}. 



2(s-2) 4 A 
Then we obtain (13.91) . For /a we set 

and the first integrand is written using the reciprocity law (13. 7p as 



^u -5/2^ (-\ - u--'"^^ (-\ + u-'/'$-^^ (- 



3 

4 V'^y du \u ) ' du^ \u 
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In-*, f i) P (1.) + !.„-3, (1.) £_, f i 



2 1 d ^ fl\\ f d'^ J\ 



\du \u) ) \dv? \u 
By similar calculation we can obtain f l3.10p . 



B Incomplete Gamma Functions 

As an example, here we only consider the integral 

/oo 
duu-^l''{d{uf -\\. 

By definition of the Jacobi theta function (11. 7p . we see that 

(ni,n2)GZ2\{{0,0)} 

/OO (yzi (yzi nQQ oo 

du u~^l^ Yl Yl e-"("'+"^^" + 4 / du u-^l^ Y. 
ni=ln2=l "^^ n=\ 



—■Kri^u 



In the first and the second integrals in the above, we set 7r(nf + n^u = w and im'^u = w, 
respectively, and replace the integral variables uhj w. Then we have 



oo oo 



us^ us^ -f /'OO '-^ -| /"OO 

L = 47r-i/2 y y^ -— -i-— / dww-'/'e-'" + A7T-'/'y- dwu,-' 

oo oo ^ oo ^ 

= 4vr-V^ Y Y 7;^^-^r(l/2, vrK + nl)) + 4.-/^ J^ ir(l/2, .n% 

ni=ln2=l '^ 1 ~'~ 2^ 



/2e-«. 



n 

71=1 



where r(z,p) is the incomplete gamma function defined by (13.151) . Combination of similar 
calculations will give the expressions (I3.16p -( 137T8|) . 
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